Affinity, Specificity, and Cooperativity of DNA Binding by Bacterial Gene Regulatory Proteins.

Nearly all of biology depends on interactions between molecules: proteins with small molecules, proteins with other proteins, nucleic acids with small molecules, and nucleic acids with proteins that regulate gene expression, our concern in this Special Issue. All those kinds of interactions, and others, constitute the vast majority of biology at the molecular level. An understanding of those interactions requires that we quantify them to learn how they interact: How strongly? With which partners? How-and how well-are different partners distinguished? This review addresses the evolution of our current understanding of the molecular origins of affinity and specificity in regulatory protein-DNA interactions, and suggests that both these properties can be modulated by cooperativity.

Crystal structures of Val58Ile tryptophan repressor in a domain-swapped array in the presence and absence of L-tryptophan.

The crystal structures of domain-swapped tryptophan repressor (TrpR) variant Val58Ile before and after soaking with the physiological ligand L-tryptophan (L-Trp) indicate that L-Trp occupies the same location in the domain-swapped form as in native dimeric TrpR and makes equivalent residue contacts. This result is unexpected because the ligand binding-site residues arise from three separate polypeptide chains in the domain-swapped form. This work represents the first published structure of a domain-swapped form of TrpR with L-Trp bound. The presented structures also show that the protein amino-terminus, whether or not it bears a disordered extension of about 20 residues, is accessible in the large solvent channels of the domain-swapped crystal form, as in the structures reported previously in this form for TrpR without N-terminal extensions. These findings inspire the exploration of L-Trp analogs and N-terminal modifications as labels to orient guest proteins that cannot otherwise be crystallized in the solvent channels of crystalline domain-swapped TrpR hosts for potential diffraction analysis.

Guest-protein incorporation into solvent channels of a protein host crystal (hostal).

Soaking small molecules into the solvent channels of protein crystals is the most common method of obtaining crystalline complexes with ligands such as substrates or inhibitors. The solvent channels of some protein crystals are large enough to allow the incorporation of macromolecules, but soaking of protein guests into protein crystals has not been reported. Such protein host crystals (here given the name hostals) incorporating guest proteins may be useful for a wide range of applications in biotechnology, for example as cargo systems or for diffraction studies analogous to the crystal sponge method. The present study takes advantage of crystals of the Escherichia coli tryptophan repressor protein (ds-TrpR) that are extensively domain-swapped and suitable for incorporating guest proteins by diffusion, as they are robust and have large solvent channels. Confocal fluorescence microscopy is used to follow the migration of cytochrome c and fluorophore-labeled calmodulin into the solvent channels of ds-TrpR crystals. The guest proteins become uniformly distributed in the crystal within weeks and enriched within the solvent channels. X-ray diffraction studies on host crystals with high concentrations of incorporated guests demonstrate that diffraction limits of ∼2.5 Šcan still be achieved. Weak electron density is observed in the solvent channels, but the guest-protein structures could not be determined by conventional crystallographic methods. Additional approaches that increase the ordering of guests in the host crystal are discussed that may support protein structure determination using the hostal system in the future. This host system may also be useful for biotechnological applications where crystallographic order of the guest is not required.

Conserved Dynamic Mechanism of Allosteric Response to L-arg in Divergent Bacterial Arginine Repressors.

Hexameric arginine repressor, ArgR, is the feedback regulator of bacterial L-arginine regulons, and sensor of L-arg that controls transcription of genes for its synthesis and catabolism. Although ArgR function, as well as its secondary, tertiary, and quaternary structures, is essentially the same in E. coli and B. subtilis, the two proteins differ significantly in sequence, including residues implicated in the response to L-arg. Molecular dynamics simulations are used here to evaluate the behavior of intact B. subtilis ArgR with and without L-arg, and are compared with prior MD results for a domain fragment of E. coli ArgR. Relative to its crystal structure, B. subtilis ArgR in absence of L-arg undergoes a large-scale rotational shift of its trimeric subassemblies that is very similar to that observed in the E. coli protein, but the residues driving rotation have distinct secondary and tertiary structural locations, and a key residue that drives rotation in E. coli is missing in B. subtilis. The similarity of trimer rotation despite different driving residues suggests that a rotational shift between trimers is integral to ArgR function. This conclusion is supported by phylogenetic analysis of distant ArgR homologs reported here that indicates at least three major groups characterized by distinct sequence motifs but predicted to undergo a common rotational transition. The dynamic consequences of L-arg binding for transcriptional activation of intact ArgR are evaluated here for the first time in two-microsecond simulations of B. subtilis ArgR. L-arg binding to intact B. subtilis ArgR causes a significant further shift in the angle of rotation between trimers that causes the N-terminal DNA-binding domains lose their interactions with the C-terminal domains, and is likely the first step toward adopting DNA-binding-competent conformations. The results aid interpretation of crystal structures of ArgR and ArgR-DNA complexes.

A residue of motif III positions the helicase domains of motor subunit HsdR in restriction-modification enzyme EcoR124I.

Type I restriction-modification enzymes differ significantly from the type II enzymes commonly used as molecular biology reagents. On hemi-methylated DNAs type I enzymes like the EcoR124I restriction-modification complex act as conventional adenine methylases at their specific target sequences, but unmethylated targets induce them to translocate thousands of base pairs through the stationary enzyme before cleaving distant sites nonspecifically. EcoR124I is a superfamily 2 DEAD-box helicase like eukaryotic double-strand DNA translocase Rad54, with two RecA-like helicase domains and seven characteristic sequence motifs that are implicated in translocation. In Rad54 a so-called extended region adjacent to motif III is involved in ATPase activity. Although the EcoR124I extended region bears sequence and structural similarities with Rad54, it does not influence ATPase or restriction activity as shown in this work, but mutagenesis of the conserved glycine residue of its motif III does alter ATPase and DNA cleavage activity. Through the lens of molecular dynamics, a full model of HsdR of EcoR124I based on available crystal structures allowed interpretation of functional effects of mutants in motif III and its extended region. The results indicate that the conserved glycine residue of motif III has a role in positioning the two helicase domains.

Multiple helical conformations of the helix-turn-helix region revealed by NOE-restrained MD simulations of tryptophan aporepressor, TrpR.

The nature of flexibility in the helix-turn-helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature-sensitive mutant showed more helical character in its helix-turn-helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild-type aporepressor are used as time-averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix-turn-helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731-740. © 2016 Wiley Periodicals, Inc.

The helical domain of the EcoR124I motor subunit participates in ATPase activity and dsDNA translocation.

Type I restriction-modification enzymes are multisubunit, multifunctional molecular machines that recognize specific DNA target sequences, and their multisubunit organization underlies their multifunctionality. EcoR124I is the archetype of Type I restriction-modification family IC and is composed of three subunit types: HsdS, HsdM, and HsdR. DNA cleavage and ATP-dependent DNA translocation activities are housed in the distinct domains of the endonuclease/motor subunit HsdR. Because the multiple functions are integrated in this large subunit of 1,038 residues, a large number of interdomain contacts might be expected. The crystal structure of EcoR124I HsdR reveals a surprisingly sparse number of contacts between helicase domain 2 and the C-terminal helical domain that is thought to be involved in assembly with HsdM. Only two potential hydrogen-bonding contacts are found in a very small contact region. In the present work, the relevance of these two potential hydrogen-bonding interactions for the multiple activities of EcoR124I is evaluated by analysing mutant enzymes using in vivo and in vitro experiments. Molecular dynamics simulations are employed to provide structural interpretation of the functional data. The results indicate that the helical C-terminal domain is involved in the DNA translocation, cleavage, and ATPase activities of HsdR, and a role in controlling those activities is suggested.

Binding and Inhibition of Spermidine Synthase from Plasmodium falciparum and Implications for In Vitro Inhibitor Testing.

The aminopropyltransferase spermidine synthase (SpdS) is a promising drug target in cancer and in protozoan diseases including malaria. Plasmodium falciparum SpdS (PfSpdS) transfers the aminopropyl group of decarboxylated S-adenosylmethionine (dcAdoMet) to putrescine or to spermidine to form spermidine or spermine, respectively. In an effort to understand why efficient inhibitors of PfSpdS have been elusive, the present study uses enzyme activity assays and isothermal titration calorimetry with verified or predicted inhibitors of PfSpdS to analyze the relationship between binding affinity as assessed by KD and inhibitory activity as assessed by IC50. The results show that some predicted inhibitors bind to the enzyme with high affinity but are poor inhibitors. Binding studies with PfSpdS substrates and products strongly support an ordered sequential mechanism in which the aminopropyl donor (dcAdoMet) site must be occupied before the aminopropyl acceptor (putrescine) site can be occupied. Analysis of the results also shows that the ordered sequential mechanism adequately accounts for the complex relationship between IC50 and KD and may explain the limited success of previous efforts at structure-based inhibitor design for PfSpdS. Based on PfSpdS active-site occupancy, we suggest a classification of ligands that can help to predict the KD-IC50 relations in future design of new inhibitors. The present findings may be relevant for other drug targets that follow an ordered sequential mechanism.

Quantum Calculations Indicate Effective Electron Transfer between FMN and Benzoquinone in a New Crystal Structure of Escherichia coli WrbA.

UNLABELLED: Quantum mechanical calculations using the Marcus equation are applied to compare the electron-transfer probability for two distinct crystal structures of the Escherichia coli protein WrbA, an FMN-dependent NAD(P)H: quinone oxidoreductase, with the bound substrate benzoquinone. The calculations indicate that the position of benzoquinone in a new structure reported here and solved at 1.33 Å resolution is more likely to be relevant for the physiological reaction of WrbA than a previously reported crystal structure in which benzoquinone is shifted by ∼5 Å. Because the true electron-acceptor substrate for WrbA is not yet known, the present results can serve to constrain computational docking attempts with potential substrates that may aid in identifying the natural substrate(s) and physiological role(s) of this enzyme. The approach used here highlights a role for quantum mechanical calculations in the interpretation of protein crystal structures.

Functional coupling of duplex translocation to DNA cleavage in a type I restriction enzyme.

Type I restriction-modification enzymes are multifunctional heteromeric complexes with DNA cleavage and ATP-dependent DNA translocation activities located on motor subunit HsdR. Functional coupling of DNA cleavage and translocation is a hallmark of the Type I restriction systems that is consistent with their proposed role in horizontal gene transfer. DNA cleavage occurs at nonspecific sites distant from the cognate recognition sequence, apparently triggered by stalled translocation. The X-ray crystal structure of the complete HsdR subunit from E. coli plasmid R124 suggested that the triggering mechanism involves interdomain contacts mediated by ATP. In the present work, in vivo and in vitro activity assays and crystal structures of three mutants of EcoR124I HsdR designed to probe this mechanism are reported. The results indicate that interdomain engagement via ATP is indeed responsible for signal transmission between the endonuclease and helicase domains of the motor subunit. A previously identified sequence motif that is shared by the RecB nucleases and some Type I endonucleases is implicated in signaling.

Three-dimensional structures of Plasmodium falciparum spermidine synthase with bound inhibitors suggest new strategies for drug design.

The enzymes of the polyamine-biosynthesis pathway have been proposed to be promising drug targets in the treatment of malaria. Spermidine synthase (SpdS; putrescine aminopropyltransferase) catalyzes the transfer of the aminopropyl moiety from decarboxylated S-adenosylmethionine to putrescine, leading to the formation of spermidine and 5'-methylthioadenosine (MTA). In this work, X-ray crystallography was used to examine ligand complexes of SpdS from the malaria parasite Plasmodium falciparum (PfSpdS). Five crystal structures were determined of PfSpdS in complex with MTA and the substrate putrescine, with MTA and spermidine, which was obtained as a result of the enzymatic reaction taking place within the crystals, with dcAdoMet and the inhibitor 4-methylaniline, with MTA and 4-aminomethylaniline, and with a compound predicted in earlier in silico screening to bind to the active site of the enzyme, benzimidazol-(2-yl)pentan-1-amine (BIPA). In contrast to the other inhibitors tested, the complex with BIPA was obtained without any ligand bound to the dcAdoMet-binding site of the enzyme. The complexes with the aniline compounds and BIPA revealed a new mode of ligand binding to PfSpdS. The observed binding mode of the ligands, and the interplay between the two substrate-binding sites and the flexible gatekeeper loop, can be used in the design of new approaches in the search for new inhibitors of SpdS.

Molecular dynamics comparison of E. coli WrbA apoprotein and holoprotein.

WrbA is a novel multimeric flavodoxin-like protein of unknown function. A recent high-resolution X-ray crystal structure of E. coli WrbA holoprotein revealed a methionine sulfoxide residue with full occupancy in the FMN-binding site, a finding that was confirmed by mass spectrometry. In an effort to evaluate whether methionine sulfoxide may have a role in WrbA function, the present analyses were undertaken using molecular dynamics simulations in combination with further mass spectrometry of the protein. Methionine sulfoxide formation upon reconstitution of purified apoWrbA with oxidized FMN is fast as judged by kinetic mass spectrometry, being complete in ∼5 h and resulting in complete conversion at the active-site methionine with minor extents of conversion at heterogeneous second sites. Analysis of methionine oxidation states during purification of holoWrbA from bacterial cells reveals that methionine is not oxidized prior to reconstitution, indicating that methionine sulfoxide is unlikely to be relevant to the function of WrbA in vivo. Although the simulation results, the first reported for WrbA, led to no hypotheses about the role of methionine sulfoxide that could be tested experimentally, they elucidated the origins of the two major differences between apo- and holoWrbA crystal structures, an alteration of inter-subunit distance and a rotational shift within the tetrameric assembly.

Interdomain communication in the endonuclease/motor subunit of type I restriction-modification enzyme EcoR124I.

Restriction-modification systems protect bacteria from foreign DNA. Type I restriction-modification enzymes are multifunctional heteromeric complexes with DNA-cleavage and ATP-dependent DNA translocation activities located on endonuclease/motor subunit HsdR. The recent structure of the first intact motor subunit of the type I restriction enzyme from plasmid EcoR124I suggested a mechanism by which stalled translocation triggers DNA cleavage via a lysine residue on the endonuclease domain that contacts ATP bound between the two helicase domains. In the present work, molecular dynamics simulations are used to explore this proposal. Molecular dynamics simulations suggest that the Lys-ATP contact alternates with a contact with a nearby loop housing the conserved QxxxY motif that had been implicated in DNA cleavage. This model is tested here using in vivo and in vitro experiments. The results indicate how local interactions are transduced to domain motions within the endonuclease/motor subunit.

Binding-competent states for L-arginine in E. coli arginine repressor apoprotein.

Arginine repressor of E. coli is a multifunctional hexameric protein that provides feedback regulation of arginine metabolism upon activation by the negatively cooperative binding of L-arginine. Interpretation of this complex system requires an understanding of the protein's conformational landscape. The ~50 kDa hexameric C-terminal domain was studied by 100 ns molecular dynamics simulations in the presence and absence of the six L-arg ligands that bind at the trimer-trimer interface. A rotational shift between trimers followed by rotational oscillation occurs in the production phase of the simulations only when L-arg is absent. Analysis of the system reveals that the degree of rotation is correlated with the number of hydrogen bonds across the trimer interface. The trajectory presents frames with one or more apparently open binding sites into which one L-arg could be docked successfully in three different instances, indicating that a binding-competent state of the system is occasionally sampled. Simulations of the resulting singly-liganded systems reveal for the first time that the binding of one L-arg results in a holoprotein-like conformational distribution.

1.2†Šresolution crystal structure of Escherichia coli WrbA holoprotein.

The Escherichia coli protein WrbA, an FMN-dependent NAD(P)H:quinone oxidoreductase, was crystallized under new conditions in the presence of FAD or the native cofactor FMN. Slow-growing deep yellow crystals formed with FAD display the tetragonal bipyramidal shape typical for WrbA and diffract to 1.2 Šresolution, the highest yet reported. Faster-growing deep yellow crystals formed with FMN display an atypical shape, but diffract to only ∼1.6 Šresolution and are not analysed further here. The 1.2 Šresolution structure detailed here revealed only FMN in the active site and no electron density that can accommodate the missing parts of FAD. The very high resolution supports the modelling of the FMN isoalloxazine with a small but distinct propeller twist, apparently the first experimental observation of this predicted conformation, which appears to be enforced by the protein through a network of hydrogen bonds. Comparison of the electron density of the twisted isoalloxazine ring with the results of QM/MM simulations is compatible with the oxidized redox state. The very high resolution also supports the unique refinement of Met10 as the sulfoxide, confirmed by mass spectrometry. Bond lengths, intramolecular distances, and the pattern of hydrogen-bond donors and acceptors suggest the cofactor may interact with Met10. Slow incorporation of FMN, which is present as a trace contaminant in stocks of FAD, into growing crystals may be responsible for the near-atomic resolution, but a direct effect of the conformation of FMN and/or Met10 sulfoxide cannot be ruled out.

Calmodulin transduces Ca2+ oscillations into differential regulation of its target proteins.

Diverse physiological processes are regulated differentially by Ca(2+) oscillations through the common regulatory hub calmodulin. The capacity of calmodulin to combine specificity with promiscuity remains to be resolved. Here we propose a mechanism based on the molecular properties of calmodulin, its two domains with separate Ca(2+) binding affinities, and target exchange rates that depend on both target identity and Ca(2+) occupancy. The binding dynamics among Ca(2+), Mg(2+), calmodulin, and its targets were modeled with mass-action differential equations based on experimentally determined protein concentrations and rate constants. The model predicts that the activation of calcineurin and nitric oxide synthase depends nonmonotonically on Ca(2+)-oscillation frequency. Preferential activation reaches a maximum at a target-specific frequency. Differential activation arises from the accumulation of inactive calmodulin-target intermediate complexes between Ca(2+) transients. Their accumulation provides the system with hysteresis and favors activation of some targets at the expense of others. The generality of this result was tested by simulating 60 000 networks with two, four, or eight targets with concentrations and rate constants from experimentally determined ranges. Most networks exhibit differential activation that increases in magnitude with the number of targets. Moreover, differential activation increases with decreasing calmodulin concentration due to competition among targets. The results rationalize calmodulin signaling in terms of the network topology and the molecular properties of calmodulin.

Raman spectroscopy adds complementary detail to the high-resolution x-ray crystal structure of photosynthetic PsbP from Spinacia oleracea.

Raman microscopy permits structural analysis of protein crystals in situ in hanging drops, allowing for comparison with Raman measurements in solution. Nevertheless, the two methods sometimes reveal subtle differences in structure that are often ascribed to the water layer surrounding the protein. The novel method of drop-coating deposition Raman spectropscopy (DCDR) exploits an intermediate phase that, although nominally "dry," has been shown to preserve protein structural features present in solution. The potential of this new approach to bridge the structural gap between proteins in solution and in crystals is explored here with extrinsic protein PsbP of photosystem II from Spinacia oleracea. In the high-resolution (1.98 Å) x-ray crystal structure of PsbP reported here, several segments of the protein chain are present but unresolved. Analysis of the three kinds of Raman spectra of PsbP suggests that most of the subtle differences can indeed be attributed to the water envelope, which is shown here to have a similar Raman intensity in glassy and crystal states. Using molecular dynamics simulations cross-validated by Raman solution data, two unresolved segments of the PsbP crystal structure were modeled as loops, and the amino terminus was inferred to contain an additional beta segment. The complete PsbP structure was compared with that of the PsbP-like protein CyanoP, which plays a more peripheral role in photosystem II function. The comparison suggests possible interaction surfaces of PsbP with higher-plant photosystem II. This work provides the first complete structural picture of this key protein, and it represents the first systematic comparison of Raman data from solution, glassy, and crystalline states of a protein.

Biphasic kinetic behavior of E. coli WrbA, an FMN-dependent NAD(P)H:quinone oxidoreductase.

The E. coli protein WrbA is an FMN-dependent NAD(P)H:quinone oxidoreductase that has been implicated in oxidative defense. Three subunits of the tetrameric enzyme contribute to each of four identical, cavernous active sites that appear to accommodate NAD(P)H or various quinones, but not simultaneously, suggesting an obligate tetramer with a ping-pong mechanism in which NAD departs before oxidized quinone binds. The present work was undertaken to evaluate these suggestions and to characterize the kinetic behavior of WrbA. Steady-state kinetics results reveal that WrbA conforms to a ping-pong mechanism with respect to the constancy of the apparent Vmax to Km ratio with substrate concentration. However, the competitive/non-competitive patterns of product inhibition, though consistent with the general class of bi-substrate reactions, do not exclude a minor contribution from additional forms of the enzyme. NMR results support the presence of additional enzyme forms. Docking and energy calculations find that electron-transfer-competent binding sites for NADH and benzoquinone present severe steric overlap, consistent with the ping-pong mechanism. Unexpectedly, plots of initial velocity as a function of either NADH or benzoquinone concentration present one or two Michaelis-Menten phases depending on the temperature at which the enzyme is held prior to assay. The effect of temperature is reversible, suggesting an intramolecular conformational process. WrbA shares these and other details of its kinetic behavior with mammalian DT-diaphorase, an FAD-dependent NAD(P)H:quinone oxidoreductase. An extensive literature review reveals several other enzymes with two-plateau kinetic plots, but in no case has a molecular explanation been elucidated. Preliminary sedimentation velocity analysis of WrbA indicates a large shift in size of the multimer with temperature, suggesting that subunit assembly coupled to substrate binding may underlie the two-plateau behavior. An additional aim of this report is to bring under wider attention the apparently widespread phenomenon of two-plateau Michaelis-Menten plots.

Environment-dependent long-range structural distortion in a temperature-sensitive point mutant.

Extensive environment-dependent rearrangement of the helix-turn-helix DNA recognition region and adjacent L-tryptophan binding pocket is reported in the crystal structure of dimeric E. coli trp aporepressor with point mutation Leu75Phe. In one of two subunits, the eight residues immediately C-terminal to the mutation are shifted forward in helical register by three positions, and the five following residues form an extrahelical loop accommodating the register shift. In contrast, the second subunit has wildtype-like conformation, as do both subunits in an isomorphous wildtype control structure. Treated together as an ensemble pair, the distorted and wildtype-like conformations of the mutant apoprotein agree more fully than either conformation alone with previously reported NOE measurements, and account more completely for its diverse biochemical and biophysical properties. The register-shifted segment Ile79-Ala80-Thr81-Ile82-Thr83 is helical in both conformations despite low helical propensity, suggesting an important structural role for the steric constraints imposed by ß-branched residues in helical conformation.

Synergy of molecular dynamics and isothermal titration calorimetry in studies of allostery.

Despite decades of intensive study, allosteric effects have eluded an intellectually satisfying integrated understanding that includes a description of the reaction coordinate in terms of species distributions of structures and free energy levels in the conformational ensemble. This chapter illustrates a way to fill this gap by interpreting thermodynamic and structural results through the lens of molecular dynamics simulation analysis to link atomic-level detail with global response. In this synergistic approach molecular dynamics forms an integral part of a feedback loop of hypothesis, experimental design, and interpretation that conforms to the scientific method.